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Real-time tumor tracking in external radiotherapy can be achieved by diagnostic (kV) X-ray imaging with a 
dynamic flat-panel detector (FPD). It is important to keep the patient dose as low as possible while main- 
taining tracking accuracy. A simulation approach would be helpful to optimize the imaging conditions. This 
study was performed to develop a computer simulation platform based on a noise property of the imaging 
system for the evaluation of tracking accuracy at any noise level. Flat-field images were obtained using a 
direct-type dynamic FPD, and noise power spectrum (NPS) analysis was performed. The relationship 
between incident quantum number and pixel value was addressed, and a conversion function was created. 
The pixel values were converted into a map of quantum number using the conversion function, and the 
map was then input into the random number generator to simulate image noise. Simulation images were 
provided at different noise levels by changing the incident quantum numbers. Subsequently, an implanted 
marker was tracked automatically and the maximum tracking errors were calculated at different noise levels. 
The results indicated that the maximum tracking error increased with decreasing incident quantum number 
in flat-field images with an implanted marker. In addition, the range of errors increased with decreasing in- 
cident quantum number. The present method could be used to determine the relationship between image 
noise and tracking accuracy. The results indicated that the simulation approach would aid in determining ex- 
posure dose conditions according to the necessary tracking accuracy. 
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INTRODUCTION 

Dynamic flat-panel detectors (FPDs) are commonly used in 
clinical practice. In extemal radiotherapy, real-time tumor 
tracking can be achieved by diagnostic (kV) X-ray imaging 
with a dynamic FPD [1-3]. There is concern regarding the re- 
lationship between image quality and accuracy of target track- 
ing, because low image quality is associated with the risk of 
increased tracking errors, although it can reduce the patient 
dose. There are a number of factors that affect image quality, 
such as X-ray tube voltage, radiation dose and patient body 
shape. In particular, reducing the radiation dose leads directly 
to an increase in image noise, and this may be one of the 
major factors reducing the accuracy of target tracking. 

Recently, several methods for measuring the temporal 
modulation transfer function (MTF) and the detective 



quantum efficiency (DQE) have been proposed, and the 
properties of FPDs used in dynamic imaging have been 
reported [4-7]. There are several reports on the perform- 
ance of electric portal imaging devices [8-9]. In a previous 
study, it was revealed that a patient dose could be reduced 
by approximately 28% by optimal settings for the low-dose 
acquisition mode with respect to image quality and dose 
[10]. However, there have been no studies regarding the 
effects of image noise on tracking accuracy. It is necessary 
to address the relationship between image noise and accur- 
acy of target tracking to keep the patient dose as low as 
possible while maintaining tracking accuracy. 

A number of researchers reported on methods for simu- 
lating reduced dose images. The technique provides 
reduced dose images by adding white Gaussian noise with 
a certain standard deviation to the original image [11-14] 
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or by adjusting noise power spectrum (NPS) and DQE 
[15-16]. In general, X-ray images have image noise due to 
statistical fluctuations in the number of incident quanta 
entering a detector {q) [17-18]; a higher quantum number 
results in less image noise. The value of q is determined as 
the reciprocal of the Wiener spectrum (WS) of the system 
(i.e. q = lAVS), which is measured in flat-field images. The 
relationship between q and pixel value in an image is also 
determined from the average pixel value in the region of 
interest (ROl), where the WS was measured. Furthermore, 
q follows a Poisson distribution. Thus, images with various 
noise levels can be simulated by changing the value of q 
and inputting the values into a Poisson random number 
generator. The simulation images allow us to evaluate the 
accuracy of target tracking at any noise level and to deter- 
mine the appropriate exposure dose. We evaluated the ac- 
curacy of target tracking in simulation images with various 
noise levels. Our purpose was to develop a computer simu- 
lation method to determine imaging conditions during 
target tracking in radiotherapy. We developed a computer 
simulation platform based on a noise property of the 
imaging system and investigated the feasibility of the simu- 
lation approach. 

MATERIALS AND METHODS 

Measurement of noise power spectrum (NPS) 
Image data set 

A set of images was generated for determination of the NPS 
of a direct-type (a-Se/TFT) FPD system (SONIALVISION 
Safire2; Shimadzu, Kyoto, Japan). The system was devel- 
oped for real-time target tracking in 3D for external radio- 
therapy, and two detectors were orthogonally-placed at 1.0 
m in the source-to-image distance (SID). In this study, one 
of the FPDs was examined. An RQA5 X-ray spectrum was 
used (HVL = 7.1 mm Al, realized with 21 mm Al additional 
filtration at 70 kV) [19]. The matrix size was 2048x2048 
pixels, the pixel size was 0.123 x 0.123 mm and the field of 
view was 25.4 x 25.4 cm. Image preprocessing consisted of 
offset and gain correction as well as compensation for 
defective or nonlinear pixels, as applied in normal clinical 
use of the detector. Pixel scaling was linear with respect to 
exposure, with a bit depth of 16 bits. 

NPS determination methods 

For determination of the NPS, three independent flat-field 
images were obtained at each of two exposure levels (six 
images in total); the exposure levels (air kerma) were 7.54 
pGy and 15.7 |jGy, respectively, for the two series. The air 
kerma values were measured free-in-air in the detector 
plane with an ionization chamber (AE-132a 2902209; 
Oyogiken Inc., Tokyo, Japan). The SID was limited to 1.0 
m in the system evaluated. The ionization chamber was 
placed 500 mm behind the detector, which was located 



approximately halfway between the X-ray tube and the de- 
tector surface. The air kerma at the detector surface was 
calculated by the inverse square distance law. 

Regions of interest (ROIs), located manually near the de- 
tector center, were 256 x 256 pixels in size, with a pixel 
sampling pitch of 0.123 mm, in the same subarea of the 
full detector area, for the two series. Average pixel values 
were measured by use of Image-J ver. 1.42 [20] in each 
ROI. The NPS was calculated according to IEC6220-1-1 
[21]. For removing long-range background trends, a 2D 
second order polynomial was fitted to each image and 
subtracted. The area of each image was divided into half- 
overlapping ROIs for each image and the results were aver- 
aged. The 2D NPS was then calculated by application of 
the fast Fourier transform to each ROI. One-dimensional 
cuts through the 2D NPS were obtained by averaging of the 
central ± 7 lines (excluding the axis) around the horizontal 
and vertical axes [22]. 

Creation of conversion function from pixel value to 
quantum number 

The q is determined as the reciprocal of the WS [mm^] of 
the system as follows [17-18]: 



In the present study, to determine q using Eq. (1), the aver- 
aged WS through all spatial frequencies in two directions 
was used as WS in Eq. (1) for each exposure level. The 
quantum number per pixel ^' was then derived as follows: 

q = q 'K ps x ps (2) 

where ps is the pixel size, which was 0.123 mm in this 
study. The average pixel value (in digital units) vs the 
number of incident quanta (in count units) was fitted with a 
linear function, y — a. + hx 

Simulation of image noise 

In this preliminary study, we simulated image noise by using 
a simplified method based on statistical fluctuations in the 
number of incident quanta entering a detector [11-14]. 
Although there is a limitation to simulating image noise with 
high accuracy, it was confirmed that the method would give 
promising results for simulating image noise under different 
dose levels. 

A tracking implanted marker (home-made metal material) 
with an acrylic plate 20 cm thick was located in clinical set- 
tings during target tracking in radiotherapy and was imaged 
at the imaging conditions used clinically in our institution, 
70 kV, 250 mA, 36 ms, and SID = 1.0 m. An averaging 
image was then created from ten images obtained at the 
above dose as a substitute for the image with vanishingly 
low image noise obtained with a large dose. 
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Pixel values were converted to quantum number accord- 
ing to the conversion function. Subsequently, the resulting 
image was weighted from 0.1 to 1.0 in increments of 0.1. 
Image noise was induced by statistical fluctuation of the 
quantum incident to the detector, which followed a Poisson 
distribution. The expected value of a Poisson-distributed 
random variable is equal to the mean and variance. Signal to 
noise ratio is defined as V<7 (= ql4q) and thus a higher 
quantum number q results in less image noise [14, 17-18]. 
Thus, to simulate image noise, the weighted images were 
input into the Poisson random number generator in each 
pixel. The output was the final resulting image with image 
noise. 

Data analysis 
Target tracking 

The targets in the simulation images were tracked by a 
template-matching technique [23]. The sum of differences 
in pixel value {R) between the search area in the next 
frame, S (x + dx, y + dy), and the template in the current 
frame, T (x, y), was expressed as follows: 

N M 

^ = ^ ^ + ^/.x, 3; + ^/y) - I 

{0<x<M,Q<y<N,-lO<dx< 10,-10 <dy< 10) 

(3) 

M and are the size of the template, and dx and dy are the 
search range. The smallest R value was obtained when 
there were more similarities in the search area and template. 
The amount of shift (dx, dy) in the search area was deter- 
mined by minimizing R, and the coordinates after move- 
ment were expressed as (x + dx, y + dy). In this study, the 
initial template was given as a region into which the target 
was inserted in the first frame. After the second frame, the 
matching region of interest in the previous frame was used 
as the new template. The size of the template was 50 x 50, 
the search range was ± 10 pixels, thus the search area was 
70 X 70 pixels, determined to cover the displacements of 
implanted targets. 

Evaluation method 

Tracking accuracy was evaluated in images at 10 different 
simulated noise levels. The implanted marker was shifted in 
known amounts by image processing, ± 3 and ± 6 pixels in 
the superior-inferior and right-left directions, respectively. 
A total of nine patterns were assessed for each noise level, 
as shown in Fig. 1. The coordinates tracked were compared 
to the known shift amounts, and the differences were calcu- 
lated as tracking errors. The maximum tracking errors were 
calculated and compared between simulation images at dif- 
ferent levels of image noise. 
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Fig. 1. Markers were shifted in nine combinations of ± 3 and ± 6 
pixels in superior-inferior and right-left directions, respectively. 

RESULTS 

NPS properties 

Figure 2 shows the WS in the horizontal and vertical direc- 
tions for two exposure levels. The results indicated that the 
present system had a stable WS through all spatial frequen- 
cies, reflecting the noise property of a direct-type of FPD. 
A high exposure level resulted in a lower WS than a low 
exposure level at all of the spatial frequencies in both hori- 
zontal and vertical directions. 



Conversion function from pixel value to quantum 
number 

The average WS through all special frequencies of two 
directions in 7.54 |aGy and 15.7 (jGy were 9.57 x 10~* mm^ 
and 5.83 x 10"^ mm^, respectively. Thus, the quantum 
numbers q at each exposure level were 1.05 x 10"^ mm~^ 
and 1.71 X 10"^ mm~^, respectively, according to Eq. (1). 
The quantum numbers per pixel c( were 1.58 x 10""* pixeF^ 
and 2.59 x 10"^ pixeF^, respectively, according to Eq. (2). 
The fitted parameters values were a = 0.0585 digital units 
and b = 541.35 digital units per count. The quantum 
number per pixel c[ was not zero even when the pixel 
value was zero due to system noise caused by the electrical 
circuit and system. 
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Fig. 2. Noise power spectra as determined for the set of flat-field images at two noise levels, (a) Horizontal 
direction, (b) Vertical direction. 
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Fig. 3. Simulation images at 10 different noise levels, (a) 
Averaging image (i.e. image without noise), (b)-(k) Images with 
simulated noise achieved by decreasing the number of incident 
quanta by 10%. Image noise increased as the incident quantum 
number decreased. 

Effects of image noise on target tracking 

Figure 3 shows the simulation images at 10 different noise 
levels. The results indicated that a lower quantum number 
resulted in more image noise. It was difficult to recognize 
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Fig. 4. Relationship between the maximum tracking error and 
ratio of incident quantum number to FPD (flat-field image). The 
average image without noise has no error, while there are tracking 
errors in the images simulated in ratio of incident quantum 
number from 1.0 to 0.1. Error bars show ±SD. (SD = standard 
deviation, n = 9). 



the location of the marker at quantum numbers <40%. 
Figure 4 shows the results regarding automatic tracking of 
implanted markers. There was no error in the average 
images without noise. The maximum tracking error 
increased with decreasing incident quantum number, as 
shown in Fig. 4. In the specific settings examined in this 
study, the tracking error tended to increase at less than half 
of the original quantum number. Error bars show the stand- 
ard deviation of nine data sets (n = 9). The range of errors 
also became larger in simulation images created with 
smaller incident quantum numbers. 

DISCUSSION 

Image noise has a big effect on visualization of an object 
with low contrast like a target in radiotherapy. The present 
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method was able to provide the relationship between image 
noise levels and tracking accuracy, although it is a specific 
parameter setting. The maximum tracking error increased 
with increases in the image noise. The range of errors also 
increased with increasing image noise. In this study, the 
tracking error gradually increased after about half of the ori- 
ginal quantum number. It was actually difficult to identify 
implanted markers in images generated by less than half of 
the original incident quantum number. Such information 
would be very useful for physicists to determine the expos- 
ure dose according to the necessary tracking accuracy. The 
present method could be applied to a different FPD system, 
requiring only the determination of the conversion function 
for that system. These results indicated the feasibility of the 
simulation approach for determination of the exposure dose 
during real-time target tracking in radiotherapy. The present 
method would be used as follows: two high-quality images 
would be taken of a patient during respiration, and then 
input to the system to simulate images with various levels 
of image noise. A target would be tracked in two images 
and the tracking error estimated for each dose level. The 
results would then allow us to determine the lowest patient 
dose possible while maintaining tracking accuracy. 

There are several limitations for this method. First of all, 
it has not been generalized enough for clinical use. More 
imaging parameters and situations need to be considered to 
provide useful information for clinical practice. In this pre- 
liminary study, however, we have demonstrated the feasibil- 
ity of the simulation approach for determination of imaging 
conditions based on the relationship between noise level 
and tracking accuracy. We believe that this simulation ap- 
proach is one possible way of optimizing imaging condi- 
tions. Second, there are the other factors that affect image 
quality, such as X-ray tube voltage, image lag, image blur- 
ring and patient body thickness. In addition, it is necessary 
to consider the noise factors apart from quantum noise, e.g. 
electrical noise and structural noise. It should be high- 
lighted in this context that the simulation approach allows 
us to evaluate the accuracy of target tracking at any noise 
level. However, all of the noise factors are not involved in 
the simulation image in the present system. One solution is 
the adaption of more accurate simulation methods [22-23]. 
Furthermore, the present system could not demonstrate rela- 
tionship between exposure dose and noise level. For clinical 
implementation, further studies are required to expand the 
system, to install a high-accuracy algorithm for noise simu- 
lation and target tracking, and to evaluate it for a real 
moving target in clinical cases. 

CONCLUSION 

The present study was performed for the development of a 
computer simulation method for determining imaging con- 
ditions during target tracking in radiotherapy. Image noise 



was simulated based on the noise property of the system, 
and the simulation was able to describe the relationship 
between image noise levels and accuracy of target tracking, 
in which the maximum tracking error increased with a de- 
crease in the incident quantum number. These results indi- 
cated that the simulation approach is one possible method 
of optimizing imaging conditions. 
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